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We demonstrate that the A'' = 1 rescattering contribution to the gluon radiation from a fast 
massless quark in eA DIS vanishes in the collinear approximation. It is shown that the n onze ro 
A'^ = 1 gluon spectrum obtained in the higher-twist approach by Guo, Wang and Zhang^^Bl jg 
a consequence of unjustified neglecting some important terms in the collinear expansion. 



1. There are several approaches to the induced gluon emission fro m fast p artons due to mul- 
tiple scattering in cold nuclear matter and hot quark-gluon plasmam^EHEl^ ^j^he most general 
approach to thi s phenomenon is the so-called light-cone path integral (LCPI) approach ^ (for 
reviews, see^^^. This formalism reproduces the predictions of the BDMPS ^ and GLV ^ ap- 
proaches in their applicability regions^S!^ (at strong Landau-Pomeranchuk-Migdal suppression 
for massless partons, and thin plasmas, respectively). However, the relation b etwe en the LCPI 
approach and the h ighe r-twist formalism by Guo, Wang and Zhang (GWZ) 12141 ^ot clear. 
The GWZ approach^Sl^l jg ^j^sed on the Feynman diagram formalism and collinear expansion. 
It includes only the = 1 rescattering and has originally been derived for the gluon emission 
from a fast quark produced i n eA DIS. The analyses ^ ^ ^ neglect the quantum nonlocality in 
production of fast partons. In^^l^l nonlocal fast quark production and gluon emission have 
been treated on even footing. However, one can show that in the applicability region of the 
GWZ formalism the quantum nonlo calit y in the quark production is not important for gluon 
emission, and the gluon spectrum of^^^^ should coincide with the = 1 gluon spectrum in the 
LCPI approach But this is not the case. The GWZ gluon spectrum predicted in"^ contains 
the logarithmically dependent nucleon gluon density, which is absent in the LCPI calculations 

m 



We will demonstrate that the approximations used in^^l ^ re allv lead to a disagreement with 
the LCPI approach^. However, contrary to the results of^^EI ^j^g correct use of the collinear 
expansion gives a zero gluon spectrum. The nonzero spectrum obtained in^SHl jg g, consequence 
of unjustified neglecting some important terms. 

2. We consider the gluon emission from a fast quark produced in eA DIS for the Bjorken variable 
xb and photon virtuality Q. The transverse momentum integrated distribution for the gq final 
state can be described in terms of the semi-inclusive nuclear hadronic tensor dW^'^ /dz (hereafter 
z = u/E, where uj is the gluon energy and E is the struck quark energy). The spin effects in the 
rescatterings of fast partons can be neglected. This ensures that the spin structure of dW^'^ /dz 
is the same as for the usual hadronic tensor Wj^'^ in eN DIS. It allows one to describe the 
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gluon emission in terms of the scalar semi-inclusive quark distribution. Neglecting the EMC 
and shadowing effects it can be written as 



dfA f r , .dfN{r) . 
— = ydrn^(r)^, (1) 

where dfN{r)/dz is the in- medium semi-inclusive quark distribution for a nucleon located at r, 
and nyi(r) is the nucleus number density. 

In the LCPI approach^^ the matrix element of the q — > gq' in-medium transition is written 
in terms of the wave functions of the initial quark and final quark and gluon in the nucleus color 
field (we omit the color factors and indices) 

{gq'\S\q) = zg J dv%, {y)-i^^ Al{y)%{y) . (2) 

Each quark wave function in (151) is written as ipiy) = ^^^^^^^=^=^^'UA0(y i yr) , where A is the 

quark helicity, u\ is the Dirac spinor operator, y^ = (y^ it y'^)/\/2. The y~ dependence of the 
transverse wave functions (j) is governed by the two-dimensional Schrodinger equation 

. d(l){y-,yT) f [{pT-9AT? + ml] +^ 

with jj, = p~ . The wave function of the emitted gluon can be represented in a similar way. The y~ 
evolution of the transverse wave functions can be described in terms of the Green's function /C for 
the Schrodinger equation ([3]). One can show that for the gauges with potential vanishing at large 
distances (say, covariant gauges, or Coulomb gauge) one can ignore the transverse potential At- 
To calculate the A^ = 1 rescattering contribution we need not use the path integral representation 
for the Green's functions (which is used at the last stage of calculations^. To obtain the A^ = l 
contribution it is enough to expand K. to the second order in the external potential. It allows 
one to describe the induced gluon emission in eA DIS in terms of the free Green's functions 
K for fast partons and the gluon correlator {A'^ {yi)A'^ {y2)) in the nucleus. Diagrammatically 
it is represented by a set of diagrams like shown in Fig. 1 in which the horizontal solid line 
corresponds to K (— >) and K* (<— ), the gluon line shows the gluon correlators, the vertical 
dashed line shows the transverse density matrices of the final quark and gluon at very large y~ . 




c d 

Figure 1. 



The typical difference in the coordinate y~ for the upper and lower ^*qq vertices (which gives 
the scale of the quantum nonlocality of the fast quark production) is given by the well known 
loffe length L/ = I/uinxb- For the nucleon quark distribution L/ is the dominating scale in the 
Collins-Soper formula^] = ^ / dy-e^^^^^+y- (iV|^(_y-/2)^+^(y-/2)|iV). For the gq final 
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state the integration over the y~ coordinate of the "y*qq vertex is affected by the integration 
over the positions of rescatterings and the q ^ gq sphtting. However, for moderate xb when 
Lj <^ Ra one can neglect the effect of rescatterings on the integration over y~ . For production 
of the final gq states with Mg^ <C the restriction on y~ from the splitting point can also be 
ignored. Indeed, the typical scale in integrating over the splitting points is given by the gluon 
formation length Lf ^ E/Mg^ which is much bigger than L/ at Mg^ <C Q^. This is valid for 
both the vacuum DGLAP and the induced gluon emission. Also, at Lj <ti Lj one can take for 
the lower limit of the integration over the splitting points for the upper and lower parts of the 
diagrams in Fig. 1 the position of the struck nucleon. Then the quark production and gluon 
emission become independent and the d/Ar(r)/(iz can be approximated by the factorized form 



dfN jr) 
dz 



TN- 



dP{r) 

dz 



(4) 



where dP/ dz is the induced gluon spectrum described by the right parts of the diagrams evalu- 
ated neglecting the quantum nonlocality of the fast quark production. 

Due to confinement the typical separation of the arguments in the gluon correlators is of the 
order of the nucleon radius, Rn- It allows one to replace the fast parton propagators between 
the gluon fields in the graphs like Fig. lb by 5 functions in impact parameter space. This 
approximation is valid for parton energy S> 1/Rn- It follows from the Schrodinger diffusion 
relation for the parton transverse motion ~ L/E. Also, the smallness of the fast parton 
diffusion radius at the longitudinal scale ^ Rn allows one to replace in other transverse Green's 
functions the y~ coordinates by the mean values of the arguments of the vector potentials in the 
gluon correlators. This approximation corresponds to a picture with rescatterings of fast partons 
on zero thickness scattering centers (nucleons). The inequality uj ^ Rn for the emitted gluon 
is equivalent to Rn <^ Lf. For this reason, in the picture of thin nucleons the contribution of the 
graphs like Fig. lc,d with gluon correlators connecting the initial quark and final quark or gluon 
can be neglected since they are suppressed by the small factor Rn/Lj. These approximations 
have been used in the orig inal formulation of the LCPI approach's (the BDMPS J- and GLV 
^ approaches use them as well). Note that (similarly to the case of the quark-gluon plasma 

each gluon correlator appears only in the form of an integral over Ay^ = V-i ~ Vi ^"^^ 
vt ~ vt ■ easily show that this ensures gauge invariance of the result (to leading order 

in Us). 

Ij^ISMI ^YiQ gluon emission in eA DIS is described by the diagrams like shown in Fig. 2. 
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Figure 2. 

The lower soft part is expressed in terms of the matrix element {A\'il:{{))A^ {yi)A^ {y2)ip{y-i)\A) , 
and the upper hard parts are calculated perturbatively. Due to conservation of the large p~ 
momenta of fast partons in the Feynman propagators only the Fourier components with p~ > Q 
are important. It means that the Feynman propagators are effectively reduce d to the retarded 
(in y~ coordinate) ones. One can show that the Feynman diagram treatment of^^Hl is equivalent 
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to that in terms of the transverse Green's functions. Indeed, using the representation 



K{yT,2,y2 \vT,i,yi ) 



dp+dpT exp [-ip+{y2 - y-i) + ipT{yT,2 - yr,i)] 



(2^)3 

one can write the retarded quark propagator as 



p-l 



2p- 



(5) 



47r 7o P 



^u\u\K{yT,2,y2 \yT,i,yi ) 

L A 



+i7 5{y^ - y^ )6{yT,2 - yT,i) 



(6) 



Here u\ and ux act on the variables with indices 2 and 1, respectively. The last term in ^ is 
the so-called contact term. It does not propagate in y~ and can be omitted in calculating the 
nuclear final-state interaction effects for fast partons. Using ^ and a similar representation for 
the gluon propagator the hard parts in the higher-twist method can be represented in terms 
of the transverse Green's functions as it is done in the LCPI treatment. We emphasize that 
the description of the hard parts in terms of the transverse Green's functions automatically 
includes all the processes in the GWZ approach (hard-soft, double-hard, and interferences in 
the terminology of^^Hi). 

3. The calculation of the diagrams like shown in Fig. 1 is simplified by noting that the free 
transverse Green's function can be written as 



K{yT,2,y2 \yT,i,yi ) = 0{y2 - yi )y](ApT(yT,2,y2 )0L(2/T,i,yi ) 



PT 



(7) 



where (j)pj,{yT,y~) is the plane wave solution to the Schrodinger equation for = with the 
transverse momentum pj-- It allows one to represent the upper and lower parts of the diagrams 
shown in Fig. 1 in the form / dy~dyT4'*i{yT,y~)4'l{yT,y~)4'qiyT,y~) where the outgoing and 
incoming wave functions have the form of the plane waves with sharp changes of the transverse 
momenta at the points of interactions with the external gluon fields. This method has previously 
been used iiJ^ for investigation of the role of the finite kinematical boundaries. All the hard 
parts evaluated with the help of the plane waves agree with that obtained in Refs.EH. 

The s um of the complete set of the diagrams contributing to the = 1 spectrum can be 
written as^Sin] 



dPir) 
dz 



j dinA{rT,n + i) 



dz 



(8) 



^•3 



Here d(j{z,^) / dz is the cross section of gluon emission from the fast quark produced at distance 
^ from the scattering nucleon. At z ^ 1 (we co nsider the soft gluon emission just to simplify 
the formulas) for massless partons it reads^ISIE] 



- ' -^"JT— H{pT,kT,z,(,) , 



dz 

H{pT,kT,z,i)-- 



J_ _ {PT - kT)pT 

Pt pSiPT - krY 




2Ez{l - z) 



(9) 



(10) 



Here the limit x — > is implicit, dG{k^, x)/dkT is the unintegrated nucleon gluon density in the 
Collins-Soper form^^, which at x <C 1 can also be written as 



dG{kl,x) 
dkx 



N 



1 



dp exp {—ikTp}V'^a{p) , 



(11) 



where o"(p) is the weh known dipole cross section. 

The coUinear expansion corresponds to replacement of the hard part H by its second order 
expansion in kx (we suppress all the arguments except for kT for clarity) 



H{kT) w H{kT = 0) + 



dH 



dk^ 



krp -\- 



(12) 



Only the second order term in (jl2|) is important, which gives to a logarithmic accuracy da{z, ^) /dz oc 
/ dpTxG{p\,x)V\^H{'pT,kT,z,^)\j:^^f^. But from (fTOl) one can easily obtain V|^-ff|^^^Q = Q 

It is also seen from averaging of the hard part over the azimuthal angle of k^ which gives 
{H{pt, kxi -z, 0) 0{kT — Pt)- Thus, contrary to the expected dominance of the region kx <Pt 
only the region kx > Pt contributes to the gluon emission, and formal use of the collinear 
expansion gives completely wrong result with zero gluon spectrum. 

The zero = 1 gluon spectrum in the collinear approximation agrees with prediction of 
the harmonic osc illat or approximation in the BDMPS^^ and approaches. The oscillator 

approximation in^l^ corresponds to the quadratic parametrization of the dipole cross section 
cr{p) = C p^. This parametrization is equivalent to the approximation of the vector potential by 
the linear expansion A~^{y~,yT + /?) ~ A^(y~,yr) + yr-^^ {y~ iVt) which can be traced to 
the collinear expansion in momentum space. The first term in the expansion in the density of 
the spectrum in the oscillator approximation in the BDMPS and LCPI approaches corresponds 
to N = 2, and the term with A^ = l rescattering is absent^ISl, Jn terms of the representation ([9]) 
absence of the = 1 contribution in the oscillator approximation is a consequence of the fact 
that in this case dG/dkx oc S{kT) (as one sees from (llip ). 

4. The vanishing = 1 spectrum is in a clear contradiction with the nonzero result ofSH. This 
discrepancy is strange enough since Eqs. ([8])- (jl0p are completely equivalent to the formulation 



ofima in the approximation of thin nucleons, and the effects beyond thi s ap proximation cannot 



be evaluated in the formalisml^HI. This puzzle has a simple solution. In^SHI the nonzero second 
derivative of the hard part comes from the graph shown in Fig. 2b (at z <^ 1). The authors use 
for the integration variable in the hard part of this graph the transverse momentum of the final 
gluon, It- The ^T-integrated hard part obtained in^ (Eq. 15 of^) reads (up to an unimportant 
factor) 

dh 



H{kT) oc 



{It - kr) 



(13) 



where 



R{y ,2/1,2/2 jT,kT) = -exp 



.y-ih - kT? - (1 - z){y^ - y^)ik^ - 2lTkT) 



1 — exp i 



(2/1 - 2/ )Q'T- kT? 
2q-z{l - z) 



2q-z{l - z) 



1 — exp 



. 2/2 (^T - kT? 
' 2q-z{\ - z) 



(14) 



is an analog of the last factor in the square brackets in (fTOj) for y" 7^ 0, yjf 7^ {y^ , yi 2 
correspo nd t o the quark interactions with the virtual photon and t-channel gluons, our z equals 
1 — z vc^^). In calculating V\^H{kT) the authors differentiate only the factor 1/(/t — kT? ■ 
However, the omitted terms from the factor R are important. After the It integration they 
almost completely cancel the contribution from the l/{lT — kT? term. Indeed, after putting y^ = 
y2 and changing the variable It — > {W + kT) the right-hand part of (fT3l) does not depend on kT at 
all. We emphasize that even without the change of the variable (fT3]l leads to V|^ff|^^^p= if 

one performs differentiating correctly which is not done in^3I31_ "J^Yiq difference between and 
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gives some nonzero contribution to the V|^i/|g_^^g suppressed by the smaU factor ^ {Rjy / Lf)'^ . 
Such contributions may be viewed as zero since they are beyond predictive accuracy of the 
approximations used in- ^ 0. 

5. In summary, we have demonstrated that the colhnear expansion fails in the case of gluon 
emission from a fast massless quark produced in eA DIS. In this approximation the = 1 
rescattering contribution to the gluon spectrum vanishes. The nonzero gluon spectrum obtained 
inEHI jg consequence of unjustified neglecting some important terms in the colhnear expansion. 
The established facts demonstrate that the GWZ approach is wrong. Its predictions for eA 
DIS and jet quenching in AA collisions do not make sense. 
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° K eeping the correction from the nonzero y 



in (|14p also does not make sense under the approximations used 
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